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Abstract. Irrotational dust solutions of Einstein's equations are suitable models 
to describe the general-relativistic aspects of the gravitational instability mechanism 
for the formation of cosmic structures. In this paper we study their state space by 
considering the local initial- value problem formulated in the covariant fluid approach. 
We consider a wide range of models, from homogeneous and isotropic to highly 
inhomogeneous irrotational dust models, showing how they constitute equilibrium 
configurations (invariant sets) of the dynamics. Moreover, we give the characterization 
of such configurations, which provides an initial-data characterization of the models 
under consideration. 
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1. Introduction 

Current theories of the formation of cosmic structure explain structure as the outcome 
of the amplification of small energy density fluctuations left after the inflationary 
epoch. The usually proposed mechanism for amplification is gravitational instability, 
which is essentially the tendency of the self-gravity of a density fluctuation to induce 
collapse. Within the framework of Newtonian Cosmology several successful techniques 
have been set up to describe this dynamical process (see [|l], and references therein). 
Some remarkable and illustrative results are: (i) The Jeans instability 0]. For small 
density fluctuations there is a critical wavelength, the Jeans wavelength Aj , such that 
fluctuations with a greater wavelength will grow with time producing collapse, whereas 
fluctuations with a smaller wavelength will be dispersed as sound waves, (ii) The Lin- 
Mestel-Shu instability 0] . Given a self-gravitating, uniform, non-rotating, and pressure- 
free distribution of matter at rest, spherically-symmetric collapse takes place through a 
series of spherically-symmetric states. However, the slightest departure from sphericity 
is systematically magnified, giving rise to an instability. In particular, if the initial 
distribution is composed of oblate spheroids, the system will approach an infinitely thin 
disk {pancake formation), whereas if it is composed of prolate spheroids, the system will 
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approach an infinitely thin cyhnder (/i/ament formation). Other similar studies can be 
found, e.g., in p. (iii) The Zel'dovich approximation This perturbative approach 
and its improvements (see, e.g., [|l], allow one to follow analytically the growth of 
density fluctuations into the non-linear regime and until caustics form. It predicts that 
the generic outcome of collapse is the formation of pancakes, although other structures 
like filaments or clumps (point-like singularities) can also form, but this is less probable. 

However, despite the fact that Newtonian Cosmology provides a wide-ranging 
picture of the dynamical process of the formation of cosmic structures, there are 
situations in which it cannot make reliable predictions. Of particular importance is the 
evolution of superhorizon inhomogeneities, where the causal structure of the Newtonian 
theory (in which the gravitational interaction is transmitted instantaneously) is not 
suitable (see, e.g., @, §|). Another interesting issue for which the relativistic approach 
is needed is the generation of gravitational waves by non-linear perturbations 0. 
Therefore, the relativistic study of the gravitational instability mechanism is well 
motivated. On the other hand, the relativistic theory involves approaches that, from 
a technical point of view, are very different from the Newtonian ones. For instance, 
in order to use numerical methods with a fluid approach we would need to consider 
an initial-value problem formulation and to build a numerical code for the evolution 
equations. In the Newtonian theory we would need to solve Poisson's equation at each 
time whereas in General Relativity we only need to solve the constraints initially and 
this ensures that they will be satisfied for all later times. 

Most of the relativistic approaches make assumptions on the geometric structure of 
the spacetime, like the imposition of symmetries (e.g., spherical symmetry), that make 
it impossible to obtain generic results like those described above within the Newtonian 
framework. In this sense, the covariant fluid approach [^, |rT| is an adequate starting 



point since it deals only with physical variables. Several studies have appeared within 
this framework. In |jl2|, Irrotational Dust Models (IDMs) were studied assuming only 
that the magnetic part of the Weyl tensor, Hab , was negligible. The dynamical analysis 
favoured the spindle as the generic outcome of the collapse [13|, in contrast to the 
Newtonian prediction, the pancake. However, the study of the dynamical consistency of 
the constraints led to a conjecture that only some known models were included [the 
Friedmann-Lemaitre- Robertson- Walker (FLRW), Bianchi I, and Szekeres dust models], 
leaving this question open. In ||15[ it was shown that if we relax the condition Hab = , 
assuming instead Hab to be transverse, D^Hab = 0, the whole set of constraints is, in 
general, unstable (a chain of infinite new constraints appears) . The main conclusion we 
can extract is that the dynamics of the gravitational instability mechanism is non-local 
(see also []16|)j and hence we must be careful when imposing conditions on Eab and 
Hab 0, |T^5 |T3 (other studies in this direction can be found in 



This suggests that the use of exact analytic approaches to understand the 
gravitational instability mechanism is quite unrealistic. Therefore, we should consider 
perturbative and/or numerical methods. The aim of this paper is to contribute to the 
development of the latter. To that end, we will consider IDMs. Thus, we assume that 
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pressure is negligible and that the matter fluid flow is irrotational. These assumptions 
are analogous to those made in the original Zel'dovich approximation Here, we 
will study equilibrium configurations of the dynamics (IDMs constitute an infinite- 
dimensional dynamical system), as we will call the invariant sets of the dynamics. We 
will see that some exact IDMs, from homogeneous and isotropic to inhomogeneous 
models, constitute equilibrium configurations and we will find their characterization. 
On one hand, this study provides a deeper knowledge of the dynamics, especially of the 
role of the gradients in the evolution equations, responsible for the non-locality of the 
dynamics, and which were neglected in most of the particular models considerer until 
now in the literature. On the other hand, the characterization of particular IDMs as 
equilibrium configurations is at the same time an initial-data characterization. Then we 
know the initial data whose development corresponds to those IDMs, which provides 
interesting examples to check numerical codes and also information about what kind of 
data we must not prescribe if we want to study new behaviours. Moreover, it can be 
used to identify attractors, repellers or asymptotic states in the dynamics. Furthermore, 
an initial-data characterization provides also a solution of the constraint equations, the 
solving of which presents in general strong difficulties, especially if we are interested in 
prescribing generic initial data. This is even more important in cosmology, where we do 
not have physically reasonable boundary conditions in order to get a well-posed elliptic 
problem from the constraints. In that respect, the information that our study provides 
can be useful for solving the constraints in other situations. 

The plan of the paper is as follows: In Sec. 0, using the covariant fluid approach and 
tetrad methods, and following |T^, we introduce a local Initial- Value Problem (IVP) 
formulation for IDMs. In Sec. ^, using this formulation, we see that a wide range of 
exact IDMs constitute equilibrium configurations of the dynamics. Moreover, we find the 
characterization of those states, which will provide an initial-data characterization of the 
irrotational dust models considered. We end with some comments and remarks in Sec. 0. 



The notation we use in this paper follows that of previous works, especially 0]. We use 
units in which 8nG = c = 1, round brackets enclosing indices denote symmetrization 
and square brackets antisymmetrization. Throughout the paper we use coordinate 
charts as well as tetrads in order to express tensorial components. The conventions 
for indices are the following: we denote spacetime coordinate indices by the lower-case 
Latin letters a, . . . ,1 = 0,1,2,3, and spacetime indices with respect to an arbitrary 
basis {eo, ei, 62, 63} by the remaining of lower-case Latin letters m, . . . ,z = 0, 1, 2, 3. 
When we choose a basis adapted to the dust velocity field, i.e., Bq = u, we will 
denote indices with respect to a spatial basis {61,62,63} by lower-case Greek letters 
a, . . . , A = 1,2,3 and spatial coordinate indices by the remaining of lower-case Greek 
letters /i, . . . , = 1, 2, 3. 
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2. An initial-value problem description of irrotational dust models 

Irrotational dust solutions of Einstein's equations are models suitable for describing the 
gravitational instability mechanism and the dynamics of the Universe in the matter- 
dominated era. Their energy-momentum distribution is completely described by the 
fluid velocity u and the energy density associated with it, p, the energy- momentum 
tensor then being 

Tab = pUaUb, u"'Ua = -l, p>0. 

The fluid flow is irrotational (^[aVbMc] = 0), i.e. u generates orthogonal spacelike 
hypersurfaces, and geodesic (w^Vfew" = 0) by virtue of the vanishing of the pressure and 
the energy-momentum conservation equations. Hence, there is locally a function r(x") 
such the fluid velocity is given by 

u = — — , u = —cLt . 

OT 

Then, r is at the same time the proper time of the matter and the label of the 

hypersurfaces orthogonal to u, S(ri) = {r(x'*) = t\ : constant} . Choosing three 

independent first integrals of n, y^{x°') {u"'day^ = 0), {t, y^} is a set of comoving 
geodesic normal coordinates. The line element in such a coordinate system has the 
following form 

rfs' = -rfr2 + V(r,y-)dy^rf2/% (1) 

where h^j_y{T,y^) are the non-zero components of the orthogonal projector to the fluid 
velocity, hab = Qab + ^a^b, whose restriction to a hypersurface S(r) coincides with its 
first fundamental form. 

In order to establish an IVP for IDMs in terms of variables with a clear physical 
and geometrical meaning, it will be very convenient to describe them in terms of 
the covariant fluid approach introduced by Ehlers |jTO| (see [0 for more details) 
and using an orthonormal basis adapted to the fluid velocity, {eo = 1^,61,62,63} 

[6m ■ Cn ^mn 

= diag(— 1, 1, 1, 1)]. We can partially fix the freedom in the choice 
of the triad {6q,} by requiring it to be parallely propagated along the fluid flow 

= M^Vbe„" = . 

This choice makes the local angular velocity vanish 

= le'^/'^ep ■ 65 = . 

The advantage of this choice is that we avoid having as a variable which is not a 
dynamical quantity since there is no evolution equation for it. 

The variables that we will use to describe IDMs can be divided into five groups 
(see p^ , ^ and references therein for details): (i) Metric variables. The components 
of the triad vectors in adapted coordinates {t, y^} [see Eq. (|I])], e^'^ . (ii) Connection 
variables. The spatial commutators, 7"/3a, defined by the commutation relations between 
the triad vectors, 

[e/3, ex] = 7"/3a4 , 7"[/3A] = 7"/3a • (2) 
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In this work we will also use the equivalent variables introduced by Schiicking, Kundt 
and Behr (see and references therein) 

7 /3A — 2a[/3(J A] + ^fi\&n <^=^ tto, — 27 a/3 ■, n — 7^^A<5 • 

(iii) Kinematical variables. The expansion (= VqW") and the shear tensor of the fluid 
worldlines. The shear is a symmetric and trace-free spatial |^ tensor 

CTab = h(^ahb)'^S/ dUc " (0/3)/la6 = ^ aUb - (0/3)/la6 . 

It is worth noting that the restriction of the quantity 

Kab={Q/^)hab + (Tab, (3) 

to a hypersurface S(r) coincides with its second fundamental form. In this work we 
will consider as variables the expansion, and the five independent components of the 
shear in a triad, (Jap , or their following representation: 

0-+ = -|o-il, Cr_ = ^((T22 - 0-33) , 0^1 = ^30-23, a2 = (Ti'i , (T3 = ^3 (T12 • (4) 

(iv) Matter variables. Apart from the fluid velocity, the only matter variable in our case 
is the energy density p (= Too). The Ricci tensor is then determined through Einstein's 
equations 

Rab - \Rgab + J^gab = PUaUb , (5) 

where A denotes the cosmological constant, (v) Weyl tensor variables. The Weyl 
tensor Cabcd describes the spacetime curvature not determined locally by matter. Its 
ten independent components can be divided into two spatial, symmetric and trace-free 
tensors 

Eab = CacbdU'^u'^ , Hab = *CacbdU'^u'^ {*Cabcd = \Vab'^^ C cdef) , 

where rjabcd is the spacetime volume 4-form. Using the analogy with electromagnetism, 
they are called the gravito-electric and -magnetic fields respectively. Whereas the 
gravito-electric field produces tidal forces having a Newtonian analogue (the trace- 
free part of the Hessian of the Newtonian potential), the gravito-magnetic field has 
no Newtonian analogue. The five independent components of the gravito-electric and 
-magnetic fields can be represented by the quantities (i?+, E^, Ei, E2, E3) and (if+, 
H_, Hi, H2, which are defined as in (^. 

The equations governing the behaviour of these quantities come from the Ricci 
identities applied to n, the second Bianchi identities, Einstein's equations (||) and the 
Gauss equations applied to the foliation {E(r)} (see |]1T| for details). In the covariant 
fluid approach the covariant derivative Va is decomposed into a time derivative along 
the fluid velocity. A"-'"},... = u'^'VcA"'"'b-, and a spatial covariant derivative tangent to 
the hypersurfaces S(r), DcA°-"'b... = /i"e ■ ■ ■ h-^bhc'^'VdA^'" f..., where A°-"'b-. is an arbitrary 
tensor field. It will be useful to introduce the spatial divergence and curl of an arbitrary 
2-index tensor Aab |^ (see for details on this notation) 

div {A)a = D^Aab , Cm\{A)ab = EcdiaD'Ab)'^ , 
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where Sabc = Vabcd.u'^ {^"'''^^def = Slh^"" eh^^ /) is the volume 3-form of the hypersurfaces 
S(r). The projection of these definitions onto a triad {Cq,} gives|| 

div(A)„ = (^5 - ^as)A\ - eJ^np^Axs , 

where da = Qa^dyt^ . Finally, we define the commutator of two spatial symmetric tensors, 
Aab and Bab, as 

[A, BU = 2A[a'B,], , [A, Bl = B]''^ = SabcA^B'^'^ . 

The set of dynamical equations can be divided into two groups: (i) Evolution 
equations, which give the rate of change of our variables along the fluid world-lines, 
(ii) Constraint equations, which are relations only between spatial derivatives of our 
variables. The form of these equations, written in an orthonormal basis adapted to the 



fluid velocity, is the following 

Evolution equations 

e„'^ = -(|e5/ + a/)e/, (6) 

iU = + ^\-) 7j]e + 26''6,,6^^iaa^^p}^l - Sap.H'^ , (7) 

e = -|02-a"V,^-ip + A, (8) 

C^a/S = — |0CrQ/3 — Cr<a^Cri3>S — Ea/S , (9) 

p=-Qp, (10) 

Eafj - cm\{H)ap = -QEap + 3a<jEpys - \p(Tap , (11) 

Hap + curl(E)„^ = -QHap + Sa^jHpys , (12) 

Constraint equations 

C^/ = ji^dsy^ - [da, dfs]y^' = , (13) 

C\ = diY{a)a-ldaQ = Q, (14) 

C'^ap = curl(a)«^ - Hap = , (15) 

C'a = diY{E)a - BapSCr^'H^ - \daP = , (16) 

C\ = diY{H)a + BapSCT^'E', = , (17) 

C^ = p-le^ + la''Paaf3-l^R + A = 0, (18) 

= EafS — ^Sal3 + <J <y P>5 — |0Cra/3 = . (19) 



X Angled brackets on indices denote the spatially projected, symmetric and tracefree part: ^(q/s) = 

- (^^A/3)5a/3 • 
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In these equations, and ^S^fs are the scalar curvature and the trace-free part of the 
Ricci tensor of the hypersurfaces S(r) respectively. They must be understood as given 
in terms of 7^^ and their derivatives, through the expression 

where F^^ = e" ■ (Vg^e^) are the Ricci rotation coefficients associated with the triad 

{e^}, related to j'^s by 5^,^^ = S^ellp + ^Pell.x + ^A.7a/3 • 

Now, we will study the evolution equations to show that they form a first- 

order symmetric hyperbolic system of partial differential equations (general accounts of 
this subject can be found in ||2^, and specializations in general relativity related with 
this work in HI, H). Following ^ we can write the evolution equations (|6|-|T2l) in 
matrix form as 



M''daU = N, (20) 
where we have grouped the variables (38 in all) into the vector U 

U iJ-J basisi U connection 1 U kinematical i U matter i U \\reyl) i l^-^) 

where the superscript stands for the transpose of a row vector, and where 

t/ftas^. = (61^62^63'^)^, 

U connection ('-^1; 02,03, nil, ni2,ni3, 7^22, ri23,'^33j , 
Ukinematical = (©5 CT-; '^2, CT-^)^ , U matter = (p) ) 

Uweyi = (£'+, El, E2, Es, H^,Hi, H2, H^)^ . 

The right-hand side of (120|) is the principal part of the system of equations, whose form 
is determined by the matrices Ai"', which in our case are (see also ||20| ) 

/ Idg ■ ■ ■ ■ \ / O9 ■ ■ ■ ■ \ 

• Idg • • ■ . O9 • • • 

M° = ■ ■ Ids ■ ■ , = ■ ■ Oe ■ • 

. 1 . . . . ■ 

V ■ ■ ■ ■ idio ; [ . . . . w'^ J 

where Hat and Oat denote the N x N identity and zero matrices respectively, and where 
the matrices are given by 



O5 A 
-A O5 



where 



O5 B 
-B O5 



V3 
2 

-1 



O5 c 
-C O5 



V3 \ 
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As we can see, A, B, and C are antisymmetric, and hence the matrices are 
symmetric. Moreover, iV is a vector which depends analytically on the variables U. 

This formulation shows that the evolution system ( pOD satisfies the conditions to 
be a first-order symmetric hyperbolic system (FOSHS) of partial differential equations 
(see [0, 0). Indeed, the matrices Ai°' are symmetric, and there is a 1-form t such 
that taAi"" is a positive definite matrix; this is immediate upon taking t = —u (for 
more details see [^). Moreover, if every 1-form t satisfying the last property is past- 



directed, the system is said to be causal. We can check that (^OD is causal (see also pO| ). 



Furthermore, as we have said before, FOSHSs admit a well-posed IVP |2^, which implies 
that given smooth initial data oU = C/|s(to) , there exists a unique smooth solution ^U. 

Now let us consider the constraints. First of all, we have to point out that not all 
of them are independent. We find the following relationships between them 

C\ = leC\ - \a^pC^P + [a, CX + div(C6), - \DX^ , 

C\ = icurl(Ci),-div(C2)„. 

The last one was already given in |3^. This means that if constraints C^a, C^, 
and C^ap are satisfied, constraints C^^ and C^q, will be satisfied automatically. Another 
important point concerning the constraints is their consistency with evolution. In 
I|, consistency was proved for analytical initial data. In [^, a general treatment 



for perfect fiuids is given and consistency has also been shown for smooth initial data, 
but the constraints have not been split according to the variables used. In our treatment, 
we can construct a closed system of evolution equations for V = {U, C = (C^))- It is 
worth noting that spatial derivatives only appear in the equations for C^a and C^a, and 
they have essentially the same principal part as Maxwell's equations. Therefore, the 
evolution equations found for V form a FOSHS. Then, given smooth initial data oU 
satisfying the constraints (C'^[oC/] = 0), and with *[/ being the only solution of (|3), it 
is clear that (*C/, 0) is a solution, and is unique. 

3. Equilibrium configurations and their characterization 

We have just seen that the dynamical equations for IDMs constitute a symmetric 
hyperbolic system, which means we can formulate a well-posed IVP for irrotational 
dust matter. As a dynamical system, one important question is the study of equilibrium 
configurations in the state space, or in a more technical terminology (see, e.g., |^), the 



study of invariant sets in the evolution. Here we will look for such configurations, 
associated with particular IDMs, finding the relations that characterize them, which 
also characterize the associated IDMs in terms of initial data. 

Let us start by introducing some ideas, concepts and the procedure that we will use 
systematically. First, the state of the system will be described by the dynamical variables 
contained in the vector U [Eq. (PT])], which we will call the state vector. Then, we will 
say that the state of the system at a time Ti is determined by U on S(ri), i.e. t/|s(Ti)- On 
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the other hand, by configuration we will mean sets of states characterized by relations 
which do not involve explicitly time or time derivatives. That is, relations of the form 

n''[U,d^U,dadfsU,...]=0. (22) 

For instance, the constraints [Equations (|13]-^)] have this form. Then a particular 
configuration whose characteristic relations ( P^ ) are preserved by the evolution will be 
called an equilibrium configuration (an invariant set). For these, 'R.^\j:(ti) = implies 
that TZ^ = for r > Ti. The constraints are again a good example. A consistent set of 
relations of the form (p2|), and defined in an open domain of the spacetime, determines 
a particular IDM in this domain. Therefore, this set of relations constitutes both a 
spacetime and an initial-data characterization of that IDM. 

An interesting example is the configuration defined by the vanishing of the gravito- 
electric and -magnetic fields: 

TZ^ab = Eab = , TZ^ab = Hab = . 

As is well-known, the vanishing of these two tensors in an open domain of the spacetime 
implies that the metric there is a FLRW metric. However, assuming that these quantities 
vanish only on an open domain of a hypersurface S(ti), Eab\ji(Ti) = -f^afels(ri) = 0, 
this does not imply the vanishing of these quantities at a different time. This is a 
consequence of the evolution equations (|lT|,|T2|): from (0) we see that when p 7^ 
the shear acts a source generating gravito-electric field Eab, and from (|T2|) we see that 
cyn\{E)ab acts as a source for the gravito-magnetic field Hab- Therefore, we conclude 
that Eab = Hab = 0, which is a spacetime characterization of the FLRW models, 
is not an equilibrium configuration. This refiects the obvious fact that a spacetime 
characterization is not in general an initial-data characterization. On the other hand, 
the vacuum configuration, defined by 

7^ = p = 0, 

is an equilibrium configuration as we can see from the evolution equation (10). 

To study whether a particular IDM determines an equilibrium configuration we will 
look for a set of relations (^2]) , R = (JZ^), satisfied by that IDM, and such that they 
are preserved by the evolution. To that end, we will consider the relations R as new 
variables and will study their evolution. For all the particular IDMs we are going to 
consider in this paper, we will look for evolution equations having the following form 

it + 0"[U, R]dJJ + P"[C/, R]d^R = X[U, R] . (23) 

and with the property 

C»"[L/,O] = 0, X[t/,0] = 0. 

As is clear, this property implies that i? = is a solution. Then, if we can show 
that this is the only possible solution, it follows that the relations in R determine an 
equilibrium configuration. One way to tackle this question is to consider the variables 
V = ([/, R) and their evolution equations (^ , |23D . In all the cases we will consider here. 
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the coefficients appearing in these equations are analytic functions of the variables V. 



Therefore, the system (^,^) for V will satisfy the conditions of the Cauchy-Kowaleski 
theorem, which does not assume hyperbolicity at all (a statement of this theorem can 
be found, e.g., in |2^), and hence, the uniqueness of the solution R = is ensured for 
analytical initial data (assuming, of course, this initial data oU is such that o-R = 0). 
In the cases where the complete system (EO,|23|) is symmetric or strongly P2| hyperbolic 



see, e.g., we will get a characterization also valid for smooth initial data. The 



constraints constitute again an example since we obtain a FOSHS for them. 



In what follows we will apply this procedure to the following IDMs: In Sec. to 



the FLRW models. In Sec. |3.2| to the spatially homogeneous IDMs and the Bianchi I 



subcase. In Sec. p.3| to the Szekeres models. Finally, in Sec. p.4| , to the IDMs in which 
the hypersurfaces S(r) are flat. In all these cases we will give the dynamical system 
for V that shows explicitly that = is a solution of the evolution equations for the 
relations (0) which determine the equilibrium configuration. However, in some cases 
this system for the variables V will not provide a well-posed IVP for smooth initial 
data. In these cases we can use the fact that the relations R are defined in terms of U, 
and hence we can get several different systems of evolution equations for R. Moreover, 
we can also use the fact that equations ( ]20|) are not coupled with the equations for R, 
and that given smooth initial data oU , there exists a unique smooth solution ^U{t, y^) 
such that *C/|s(t-„) = oU . Then given a system of equations for R, we can substitute U 
for in them so that we get a linear system of partial differential equations 

R+Q^[Ty]d^R = Y[Ty,R], 

where all the coefficients are smooth. If we can obtain a system like this such that 
it admits the formulation of a well-posed IVP for smooth initial data, and if we can 
show that = is a solution, then we would have shown that it is the solution. 
Therefore, our characterization would be also valid for smooth initial data. We have 
applied this procedure successfully to some of the cases enumerated above. For the 
sake of brevity, we will not write down explicitly these alternative systems, especially 
since in most cases they can be obtained by making a few changes in the systems (^) 
for R. Finally, it is important to remark that the initial-data oU in each case must 
satisfy, apart from the relations (0), the constraints. Taking this into account, the 
characterizations of the different equilibrium points that we are going to give will consist 
of the minimum conditions that, added to the constraints, imply the relations defining 
such configurations. 

3.1. FLRW dust models 

The FLRW models are the standard cosmological models. They are motivated 
by the so-called Cosmological Principle in the sense that they are homogeneous and 
isotropic cosmological models (they have a six-dimensional group of motions). In the 
case of dust, the line element can be expressed in terms of elementary functions as 
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follows (see, e.g., p^ ): 



ds^ = -dr^ + 




(l + fr^) 



where 



R = (E/3)(cosht- 1), 
R = t\ 

R = (-E/3)(l-cost), 



r = 



r = 



r = 



{E/3){smht -t) , 
{-E/3){t - sint) 



if A; < , 
if = , 
if A; > , 



E is an arbitrary constant such that sign(£') = — sign(/c), and the hypersurfaces 
S(r) have constant curvature {^R = QkR^^). There are several covariant spacetime 
characterizations of the FLRW models, and in particular of the dust subcase (see, 
e.g., [^). Here, we will be concerned with the following two characterizations: 

(i) The kinematical characterization. In terms of the kinematical quantities a dust 
FLRW model is characterized simply by the vanishing of the shear 



(ii) The Weyl tensor characterization. The FLRW dust models can be characterized 
by the vanishing of the Weyl tensor, or equivalently, by the vanishing of the gravito- 
electric and -magnetic fields 



From these spacetime characterizations we can see that the FLRW dust models 
constitute an equilibrium configuration in the dynamics of the IDMs, finding at the 
same time an initial data characterization for these models. To that end, we need to 
consider the following relations 



As is clear, = is a solution, and since the equations for {U, R) constitute a 
FOSHS, it is the only one. To sum up, we can say that if the initial state belongs 
to the configuration characterized by the vanishing of the shear and the gravito- electric 
field, the system will remain in this configuration. Moreover, the resulting spacetime 
belongs to the FLRW class. In this statement we have not taken into account the 
vanishing of the gravito-magnetic field because, as we have said before, we assume that 
the constraints hold. Apart from the vanishing of Hab, the constraints imply 



Eab — Hab 



0. 



ab = O'ab — , IZ ah = Eab — , TZ ab = Hab — . 

Their evolution follows directly from (|6|-p^) 




D,e = 0, DaP = Q, 

that is, 6 and p must be constant on the hypersurfaces S(r). 



(24) 
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We can find another initial-data characterization of the FLRW models, equivalent to 
that given above, in terms of geometrical quantities only. It is based on the equations (|^) 
and (|T8|JT9|) , and the relations that characterize it are 

n\b = - \KbK = , n\b = ^Sab = , n\b = curKi^),^ = o . 

As before, TZ^ab, Tl'^ab and the constraints imply T^^ab- The remaining conditions imposed 
by the constraints, which are equivalent to (p4|), are 

K\t.{t) = constant , ^-R|s(r) = constant . 

All these relations together, tell us that the initial hypersurface S(ro) must be of constant 
curvature and the second fundamental form (which describes how it is immersed in 
the spacetime manifold) must be proportional to the first fundamental form, with the 
proportionality factor the trace K of Kab , which is a constant. 



3.2. Spatially homogeneous IDMs 

This well-known class of spacetimes is characterized by the existence of a three- 
dimensional group of motions {G3) simply transitive on the hypersurfaces orthogonal 



to the fiuid velocity, {S(r)}. These models were studied systematically in ||2^, where a 
classification was given. Following this work we can see that a spacetime characterization 
of these models is the following: "a perfect-fluid cosmological model contains a three- 
dimensional group of motions simply transitive on three-surfaces orthogonal to the fluid 
velocity u if and only if there exits an orthonormal basis {eo = u, Cq,} such that the 
commutators satisfy 

The triad {Gq,} spans the surfaces of transitivity of the group. 

Starting from the above characterization we will prove the following statement: // 
the initial state is such that the quantities oj0s> o©; o-nd oCTa/s are constant, the system will 
stay in this configuration. The cosmological models constructed from the development 
of initial data satisfying these conditions will correspond to spatially homogeneous 
(Bianchi) dust models. To show this, let us consider the following relations 

V = ^a7^<5 = , n\ = d^Q = Q, n\^s = d^aps = , (25) 

TZ'^a = 9aP = , TZ^al35 = daEjSS = , TZ^ af35 = = . (26) 



Whereas the relations in ( pq ) will implement the characterization, the relations in (Pq) 
are auxiliary conditions needed for the proof, but they can be derived from (^) and the 
constraints (13-0). Using the evolution and constraint equations and the commutators 



of spatial derivatives (0) we obtain the following system of evolution equations for R 

n^\p5 = -(^©5a^ + CT^')n'\f55 + (|©5^[/3 + [f5)n' ^ \a\5]e + 7,^[^(H'^]7^'a + ^%|5]') 
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+ 3(7(^^7^^ + 'la\£^'^{l3Tl^\e\5)K + ^AK(/3^^^^|a|<5>'' , 
^^a/3<5 = —HpslZ^a + 3-f/'(/3^7^^|Q,|5)e — (fO^a^ + el35 + 3cr(/3^7^^|a|5>e 

As we can see, only the two last equations contain spatial derivatives, and their structure 
is essentially the same as that of equations ([TT| , p!2D . Then the extended system for (t/, R) 



is a FOSHS, with i? = the only solution. This shows that Bianchi models constitute 
an equilibrium configuration characterized by ([25|) . 

This result provides an initial-data characterization of the dust Bianchi models 
within the class of IDMs, but it says nothing about the particular classes of Bianchi 
models (see for a classification). We can study which kind of equilibrium 



configurations particular Bianchi models may constitute by refining the characterization 
given in expressions (^,^). To that end, we will have to add more conditions depending 
on the class of Bianchi models under consideration. As an example, let us consider the 
Bianchi I dust models in which the G3 group is Abelian and the line-element is 



given by 

3 

ds' = -dr^ + ^ r2p°(r - ,r)2(2/3-p.)(^^a)2 ^ ^27) 

where qT, are constants such that Yl\=i'Pa — Zla=ii^a = 1- These models can be 
characterized by adding the following conditions: 

IZ^ a = eta = , TZ^aP = ^a/S = , TZ^ a/S = Ha/3 = . 

Here TZ^a/s is a consequence of TZ^aps [Eq. (P^D], TZ'^ap, TZ^ap and the constraint ([TB|). 
We can obtain the following evolution equations for these quantities: 

iz\ = -\en\ + - |£/V/7^8^A , 

n\p = - a<„^7^8^>A + 6apa^'n\5 + \aapn^\ + aA<"£^>^'7^^6 - n\p , 

iZ^ a(3 = —'dTZ^ afS + 3(T<a'^7?.^/3>A — <aT^^ \X\I3>5 ■ 

The whole system of evolution equations for the quantities U and 7Z^ (/ = 1, . . . , 9) still 
constitutes a FOSHS, and 7Z^ = is a solution, the only possible one. Therefore, we can 
say that conditions (2t) and {TZ'^ aPiTZ^ ap) determine an equilibrium configuration in the 



dynamics of the IDMs, which corresponds to the Bianchi I dust models. In a similar way 
and using the information given in |^ we can study which equilibrium configurations 



come from the other classes of Bianchi models. 
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The characterizations given above do not have a covariant translation. However, in 
some cases it is possible to get a covariant characterization, an example being Bianchi 
I IDMs. As in previous cases we start from covariant spacetime characterizations. For 
Bianchi I models two such characterizations follow from |3^. The first one is determined 



by the vanishing of the spatial covariant derivative of the shear 

The second one says that an IDM belongs to the class of the Bianchi I dust models if 
and only if the hypersurfaces S(r) are flat 

^Rab = , 

and Eab (or equivalently aab) is non-degenerate (algebraically general). Taking this into 
account, a covariant initial-data characterization of the Bianchi I dust models is given 
by 

n\bc = DaKbc = Daiabc + iQhc) = , n\b = ^Rab = . (28) 
The evolution equations for these relations form a closed system: 

Ti-^abc = —KlZ^abc — KhcTZ^a/ — K/TZ^dbc — "^Ki/^IZ^ c)ad + '^Kd{bT^^'^ c)a 

n\, = -Kn\, + K^a'n\y - \Kabn\- + omllecaian^"] 

+Kb [\Kn\' - K'^n^d) , 

where Kab must be understood as given in terms of the expansion and shear (^). Again, 
-R = is a solution, but now the system of evolution equations we have for (t/, R) is 
not a FOSHS. However, using the procedure described above we can show that R = 
is the only possible solution, showing that (^) constitutes a stationary configuration 
associated with the Bianchi I models (| 



3.3. Szekeres dust models 

The Szekeres dust cosmological models [^| are inhomogeneous models. As was shown 



m 



39| , they do not have in general any Killing vector field. From this point of view they 
can be considered to be generic cosmological models, but from other points of view, as 
for instance the algebraic classification of spacetimes, they are not so generic since they 
are Petrov type D models, a special algebraic case. The line element can be written in 
the following diagonal form 

= -dr' + J:~'U{T,y'){dy'r + A\T,y'^) [{dy^ + k' {y' , y') {dy'f] , (29) 

where S coincides with the repeated shear eigenvalue cr^ = a^. From their spacetime 



characterization (see |^ and references therein) 



[a,E]a = 0, Hab = , Eab (or (Tab) degenerated , (30) 
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"^^-,1,2,3 = 




,^2,^3) = 0, 




T^^af3 = Ha/3 = 





"^^+,-,1,2,3 = {H^, H^,Hi, H2, H3) = 


= 0, 


7l^al3 = Map = 


= 


^ 7^V,l,2,3 = (M+,M_,Ml,M2,M3: 


= 


= Nap = 





^ 7^^,_,l,2,3 = (iV+, iV_, Ni, N2, Ns) = 


= 0, 


have introduced 


the 


following definitions Mab = cuT\{E)ab , 


and 



we can see that the evolution equations (§-|T^) become local equations [0, that is, 
the spatial derivatives disappear (there are no fields propagating in these spacetimes). 
Hence, the evolution equations are ordinary differential equations. 

We can obtain an initial-data characterization of these models starting from the 
spacetime characterization (pOl). Taking into account the form of the third condition 
in (POP, the use of an orthonormal basis adapted to the fluid velocity simplifies 
considerably the relations we will use to get the characterization. These relations are 

■7^^-, 1,2,3 = (CT-, C^l, Cr2, era) = , 7^\l 22,33,23 = i^U^ ^22, ?^33, ^23) = , (31) 

(32) 
(33) 
(34) 
(35) 

; definitions Mab = cuT\{E)ab , and Nab = 
cuT\{H)ab, which become new constraints 

C\b = cni\{E)ab -Mab = 0, (36) 

C\b = cm\{H)ab -Nab = 0. 

For these quantities we did not give evolution equations, but we can find them using 
the evolution equations (|6|-P^. The result can be expressed as follows: 

Mab = -I^Mab + |fT<a1fX, H]byc ' \pHab + 3E^a'Hb>c + Cm\{N)ab + T\b , 
Nab = -iQNab + SH^a'Hb>c " curl(M),b + J^\b , 

where 

J^^ab = ^£cd<aEb>^Dea'^^ + ^ecd<aO'b>^DeE'^^ + 3cnT\ {a-E)ab — Cre^£cd<aD^Eb>^ + CCab, (37) 
^^afe = l£cd<aHb>''DeCr'^'' + 3cUrl (a ■ H)ab — CTe^ecdKaD'^Hb^''' + CC ab , 

and 

(cr ■ E)ab = Cr^a^Ebyc, {ci ' H) ab = <a Hb>c ■ (38) 

Moreover, the quantities CCab and CC' ab represent combinations of the constraints, which 
vanish when the constraints hold. Here, we have used the well-known fact that adding 
combinations of the constraints to the evolution equations does not change the dynamics, 
but they can be useful for technical reasons and can modify the type of differential 
equations that we obtain. We will use this fact to find the initial-data characterization 
of the Szekeres models. Now we are ready to compute the evolution equations for the 
quantities TV . They can be expressed in the following form 

ie_ = _|(e + a+)7^l_ - ^(^27^^2 + (j^n^) - n^. , 
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7^^ = -|(e + a+)7^^ - j^^cx^n^ - n\ , 

n\ = -|(20 + (T+ + V3cx_)n\ - ^^(ri'JZ^ - n\ , 

-R}; = -|(20 - a+ + V^a^)n\ - ^al7^l2 - n% , 

Tt'ii = -|(e + a+)7^^l + ^^T_(7^222 - 7^233) 

+ ^{n^2n\ + nl37^l2 + 2r^237^^) + ^(a37^S - a27^l2) + , (39) 

7^'22 = 1(26 - a+ - V^Or_)7^222 - ^(v^^+ " ^-)(7^^1 - 7^233) 

+ ;^(^l27^'3 + n237^^ + 2r^l37^l2) + -^ia,n\ - a37^l3) - |(7^^ + V37^4-) , (40) 

= -1(26 - (T+ + v^^T_)7^233 - ^(v^a+ + (T_)(7^^l - n\2) 
+ 75(^l37^^2 + ^237^^ + 2nun^3) + 75(a27^l2 - ain\) - \{n\ - Vsn^^) , (4i) 

7^'23 = -|(6 + a+)7^223 - i73(^l37^S + ni2n\ - ni{R\) 

-^(2al7^l_ - a2n\ + a37^l2) - ;757^^ , (42) 

n^- = -(6 - a+)7^=^_ - ^(a27^32 - a^n\) + (e+ - ip)7^l_ + n^. , 
n\ = -(6 - a+)n\ + ^((727^33 + a37^32) + {e+ - \p)n\ + n\ , 

n\ = _[e + i(a+ + V3(T_)]7^32 + ^(a37^^ + al7^33) 

- i(E+ + VsE_ + p)n\ + 7e62 , 

7^'3 = -[6 + |(a+ - V^a^)]n\ + ^(a27^^ + al7^32) 

-\{E+-V^E^ + p)n\ + n%, 

7t^/3 = -|67^^/3 + |a<„'5[a,7^4]/3>5 - + 3£;<,'^7^4;3>^ + curl(7^6)„^ + J^\p , 

7^'a;3 = -|67^6,^ + m^Jn^pys - |div(a)'5£^^<„7^4^>^ + 3curl((T • 7^4)„^ 
-a/£^,<„D^7^^^>'^ - curl(7^^),^ . 

The main point is to show that -R = is a solution for this system of equations. As 
we can see, this question depends only on the form of the different terms that make 
up the tensor J^^ap [see Eq. (|37|)]. The explicit expressions for its components are 
given in Appendix A , where the combination of constraints CCab has been chosen in the 
following way: CC+ = CC„ = CCi = 0, and [see Eqs. (|T5|) and {^^] 

CC2 = \E+C\ + ^a+C^2 , (43) 

CC3 = |^+C23 + ^a+C73. (44) 
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With this choice we have shown in [Appendix A| that J^^a/s = when TZ^ = 0. In the 



case of the evolution equation for 71^ a/3 we have not added any constraint, i.e., we have 
chosen CC'ai3 = 0. To sum up, we can find a system such that 71^ = is a solution. The 
whole system of evolution equations for ([/, R) is not a FOSHS system, but following the 
discussion above in this section we can modify the system for R to show that TZ^ = 
is the only possible solution. The conclusion is given in the following statement: // 
the initial state of the matter fluid belongs to the configuration defined by the following 
covariant conditions 

[cr, E]a = 0, Hab = 0, Eab {or aab) degenerate, cni\{E)ab = , (45) 

then it will belong to the configuration for all later time. It is worth noting that 
conditions (|3ll - ^ ) have a covariant translation. Apart from defining an equilibrium 
configuration, conditions (^51) provide an initial-data characterization of the Szekeres 
models (the Szafron models 0] with constant pressure in the case of a non-zero 



cosmological constant). The remaining conditions are auxiliary, since they follow from 
(^) and the constraints, or they can be obtained just by fixing the remaining freedom 
in the choice of the triad {e^} ■ 

On the other hand, we can have a covariant characterization equivalent to (^) 
expressed only in terms of geometrical quantities. It is given by the following relations 

[K,^R]a = 0, curl(ir),, = 0, Cab = 0, (46) 

where Cab is proportional to the conformally-invariant Cotton- York curvature tensor 
associated with the hypersurfaces S(r), The proportionality factor is given by 

^ab ^ 1^5/3 fjab ^ ^ ^ [^^^ 1/2 ^ 

One can check that Cab can be expressed in terms of the trace-free part of the Ricci 
tensor of the hypersurfaces S(r) as follows 

Cab = CUrl{^S)ab ■ 

Cab has essentially the same properties as the Cotton- York tensor: i) It is spatial, 
symmetric and trace-free, ii) It is transverse (divergence-free): DaC""^ = 0. iii) It 
vanishes if and only if the hypersurfaces S(r) are locally conformally-flat. The vanishing 



of this tensor has been sometimes related with the absence of gravitational waves ||42 
For the Szekeres models this interpretation seems to be correct since Cab vanishes, and 
as we have said before, there is an absence of propagation signals in these models (the 
evolution equations are completely local). 

The first two conditions in (^) are clearly a direct translation of the first two 
conditions in (^5]). Therefore, to show this equivalence we only have to show that the 



third conditions are equivalent. To that end, we can use the constraint (^) to write 
Cab in the following form 

Cab = CUTl{E)ab " |0i^afe + (T'^ eScdKaD^'crby'' + lecd<a(^b>'' D ^Cj'^'' 

= cni\{E)ab-\QHab + rab. (47) 
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Then the problem reduces to checking that the terms in Tab vanish when the relations 



hold. From (P7|), Tab contains only the shear and its first derivatives and hence 
we only have to check that its components do not contain any term of the type a+da(J+ 

la , , or ni3 . Since the terms contained in Tab have the same 



or r (T^ , where F 



structure as those in J-'^ab, we can use the calculations done for J-'^at (see [Appendix 



For the components T^ , T_ , and 7^ , it is a matter of algebra to check that the 
terms of the type mentioned above do not appear, and for the components 7^ and T^ 
we need to use the components 2 and 3 of the constraint C'^ab [Eq- dH)] to see that 



they cancel. These calculations finish the proof that the conditions (^61) provide an 
initial-data characterization of the Szekeres models equivalent to that given in (|45|) . 



3.4- ID Ms with flat hypersurfaces 

The spacetime characterization of these models is obvious; they are defined by the 
vanishing of the Ricci tensor of the hypersurfaces S(r), i.e., 

^Rab = 0. (48) 

From a geometrical point of view these models can be considered as inhomogeneous 
generalizations of the standard Einstein-de Sitter cosmological model [^. Their 



relevance is enhanced by recent observations [0], which suggest that the spatial 
geometry of the local universe is close to flat. For the case of a vanishing cosmological 
constant (A = 0), the IDMs satisfying (|48|) were determined in There are two 

classes of solutions, both having a vanishing gravito-magnetic field, Hab = 0. The first 
one consists of the Bianchi I dust models (^Tf ), which are Petrov type I models. The 
second class, which has a degenerate gravito-electric tensor (Petrov type D), is the 

Since we 



subclass of the Szekeres models [3S] satisfying 



which were given in [45 



have already studied the first class, we will focus here on the second class, which is 
composed of two families of solutions. The line element of the first family is given by 

ds^ = -dr^ + [(1 + Ay + Bz)t + Cf r-^dx^ + A {dy^ + dz^) , (49) 
where A, B and C are arbitrary functions of x. The line element of the second family is 



ds^ = -dr' + 



(^-or)3 



2 d^ioT) 

3 5,(ln?7) 



dx^ + U^{t - oT)^{dy^ + dz^), (50) 



where 



u = vw . 



w 



2 I 2 

y +z 



+ 2hy + 2cz + d} 



and (a, 6, c, d), V and oT are any functions of x such that ad — h^ — c^ = 1. Remarkably, 
despite the particular character of these two solutions, both share the property of having 



no Killing vectors fields in general, like the generic Szekeres metric (p9[). 



Equilibrium Configurations in the Dynamics of Irrotational Dust Matter 



19 



Now we will show that the IDMs ( ^ , pOD constitute an equilibrium configuration in 
the dynamics of IDMs, and we will also find its characterization. To that end, we need 
to consider the following relations: (|3lD, (|33D and the alternative to (0), 

n\b = ^Rab = ^ = 7^3 = = and 7^^„,> = S\b = 'Sab = o . 

The evolution equations for them are given by the equations ( PU| - ^ ), which do not 
change, and by the equations 

aai3 = — |0'53q,^ + a^a^S^i3^s — ajs'^^ + curl(7^^)Q,^ , 

where T^js is defined by Eq. (^T]). Therefore, we have a closed system for which TZ^ = 
is a solution, which can be shown to be unique. Then we can say that if initially the 
spatial geometry is flat, the gravito -magnetic field vanishes and the shear is degenerate, 
the system will stay in this configuration. Again, we have been able to extract a covariant 
initial-data characterization from the relations that define the equilibrium configuration. 

4. Remarks and discussion 

In this paper we have considered the general dynamics of IDMs from an Initial- Value 
Problem perspective. The equations governing these models constitute an infinite- 
dimensional dynamical system due to the presence of spatial gradients. We have 
seen that a variety of well-known models, from the homogeneous and isotropic FLRW 
spacetimes to the inhomogeneous Szekeres cosmological models, constitute equilibrium 
configurations of the dynamics (invariant sets). It remains to find the characterizations 
of some important IDMs, as the remaining Bianchi classes of IDMs and the Lemaitre- 



Tolman-Bondi dust models (the initial-value problem approach made in |^ could lead 
to such characterization). Moreover, we have found the relations characterizing these 
configurations, which in most cases can be expressed in a covariant way. At the same 
time, these relations provide an initial-data characterization of the IDMs considered. 
The information obtained together with the study of the evolution of the different 
physical quantities clarifies what the dynamical content of these IDMs is, as well as 
the role of the different variables and the spatial gradients. 

This work must be considered as a first step towards the implementation of a 
more general programme to understand the full dynamics of the gravitational instability 
mechanism. In that respect, two interesting issues to be elucidated are the question of 
how a triaxial configuration will evolve in collapse (among the models treated in this 
paper, only homogeneous models have an algebraically general shear and gravito-electric 
field) , and the role of non-local effects due to the appearance of the curls of the gravito- 
electric and -magnetic fields. Such a programme would entail a method to prescribe 
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generic initial data. This means finding solutions of the constraints for any possible 
initial distribution of matter. The programme would also include a numerical code to 
evolve this initial data, i.e. to construct its development, and capable of extracting the 
physical information. The development of these two points, which does not seem an 
easy task, could be complemented by the information coming from studies based on 
perturbative methods (perturbations of a background or iterative approaches such as 
the long wave-length approximation scheme). 

The results of this paper provide useful information for the development of such 
a programme. First, they can help us in the prescription of initial data: we know 
the role of the variables and what kind of initial data we have to avoid in order to 
study new behaviours. Moreover, since the characterizations given are solutions to the 
constraints, this information can be helpful in solving them for other different cases. 
On the other hand, we can use the characterizations to check numerical codes and to 
identify attractors, repellers and asymptotic states of the evolution. 

Finally, it is worth noting that the present work can be extended to models with a 
more general energy- momentum content. In this sense, a straightforward generalization 
would be the study of irrotational perfect-fluid models, including in this way the effects of 
the pressure and hence, the inclusion of sound propagation. Furthermore, it would also 
be interesting to consider fluids with rotation to study their effects in the gravitational 
instability mechanism. This would require introducing a spacelike foliation of the 
spacetime not related with the fluid velocity, since in that case it does not generate 
orthogonal hypersurfaces. 
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Appendix A. Explicit expressions for the components of J-'^ab 

To complete the proof of the characterization of the Szekeres cosmological models in 
terms of initial data given in subsection we need to show that the quantity J-'^ab 
[see Eq. (^7])] vanishes when TZ^ = . To that end we have to work out the components 
of J-'^ai3, the projection of J^^ab with respect to the orthonormal triad {Ca}- The terms 
CCab and CC'ab, denoting combinations of the constraints, have been chosen as follows: 
CC'ab = , CC+ = CC_ = CCi = 0, and CC2 and CC3 are given in Eqs. (11,11 . 

In order to simplify the large expressions for the components of J-'^ab we will use 



Equilibrium Configurations in the Dynamics of Irrotational Dust Matter 21 

the following notation: For any components and A_ of a tensor A^/s 

A± = V3A+±A_, A± = {A+±VSA^)/V3. 
The components of the tensor (cr ■ E)ab = Vab [see Eq. (pSl)] can be written as 
V+ = -\{a+E+ - (T_7^3_) - lia^nS + a^Tl^ - 2a^n\) = -\a+E+ + Q+ , 

Pi = |(a+7^^ + E+n\) + ^((737^32 + ^271^3) , 

As we can see, the quantities Q+ , , Pi , "^2 , and P3 vanish when 71^ j = IZ^ i = . 

Taking all these definitions into account, the components of T^ab can be written as 
follows 

•^V = i {[(^2 - 302)5-+ + ni30-_ + (nil - ^^33)0-2 + nuCTi - n23C^3]7^^2- 

[(^3 - 3a3)cr_ - ^12(3-+ - (nil - ?^22)o-3 " ^^130-1 + ?^230'2]P-^3 + 

2[aia[3 + d[s(J\i\ - 3(aia[3 + ala[3)]7^=^2] + {(t ^ E ^ TZ')} + 
^ {(^3 - a3 + ni2)p3 - (^2 - a2 - ^13)^2 - 2^23^1 - 3P+7^^l- 

\ {a^\d'iV}2 + (ai + r^23)7^^ + (03 - r^l2)7^33 + |^_(7^^l - 7^222 + 
d_[dsn\ + (ai - n23)7^^ + (as + ni^)n\ - i^+(7^^l + 7^222 - 7^233)] - 
[2(9(2^3] + V3n23^+ + aiE_) + (02 + ^13)^3 - (03 - ^12)^2 + - 
[9i^3 - (a2 - ni3)^+ - (03 + ni2)Ei + - ^22 - ^33)-E2]P.^2 + 

[d-s_E2 - (03 + ni2)E_ - (02 - ^13)^1 - K^ii - ^22 - n33)^3]P-^3} , 

= - ^ {[(^2 - 3a2)cr+ + r2i3(T_ + {Uu - 7233)0-2 + ^120-1 - n2'i(y'i]Tl^ 2 + 
[(^3 - 3a3)cr_ - ni20-+ - (nil - ?^22)cr3 " ^^130-1 + ?^230-2]P-^3 " 

4[(5(2 - 3a(2)(T3) - ^(^1 - 3ai)cr+ + n23cr_ - |(n22 - n33)(Ti + 
+ 1(^130^3 - ^^120-2)]^-^ + 2[(ai - 3ai)o-(2 + (9(2 - 3a(2)o-|l|]7^33) + 
(o- ^ E,TZ^ ^TZ^)] - 

I {(^3 - a3 - 3ni2)P3 + (^2 - 02 + 3ni3)P2 - 2(9i - ai)Pi- 

v^p+(7^'22 - n\,;) + P-(7^^l - 27^222 - 2n\^)} - 

^ {^(T+[9l7^^ + (a2 - nl3)7^32 + (03 + r^l2)7^33 - (nn - ^22 - r^33)7^='-] + 
a+[927^'2 + (ai + n23)7^^ + (03 - nl2)7^=^3 + |^-(7^^l - n^22 + 7^233)] + 
d.ld^-R?^ + (ai - n23)7^^ + (02 + nl3)7^=^2 - |^+(7^^l + n^22 - 7^'33)] + 
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[2(9(2-53) + V^aiE+ + ^23-5-) - (^22 - niz)Ei - (03 - ni2)-B2 - 
(02 + niz)E^\R}^ + [(ai + 2ai - 27123)^^3 - (293 + 03 + n^2)Ei - 
(02 - ni3)-B+ + 4(a2 + ni3)-B_ - |(nii + 8/122 - r^33)-£^2]7^S + 
[(^1 + 2ai + 2n23)£^2 - (292 + 02 - ni^)Ei - (03 + 7112)^- - 4(03 - ni2)E^ + 
\{nii - 7122 + 3r^33)£^3]7^S} , (A.l) 

•^^1 = ^ {[(^2 - 302)5-+ + ni3a-_ + (nil - ri33)c^2 + rii20-i - n230-3]'^^3- 
[(^3 - 3a3)a-_ - 71120-+ - (nil - n-22)crz - ni3(7i + n-22,(T2\R? 2 - 
4[(a(2 - 3a(2)(73) - :^(9i - 3ai)(7+ + 7^230•_ - §(7122 - n2,^)ai + 
\{ni^(Tz - 7^l2(72)]7^^_ - 2[(ai - 3ai)(7[2 + (9[2 - 3a[2)(7|l|]7^^3] + 

{(T ^ E ^ n^)] + 

I {(93 - 03 - 3ni2)7'2 - (92 - a2 + 3ni3)7'3 + 2(9i - ai)V-- 
2^/3P+7^223 + (nil - 2n22 - 2n33)Pi} + 

^ {^(T+[9l7^'^_ + (02 - nl3)7^33 - (aa + nl2)7^■^2 + (nn - 7222 - n33)7^^] + 
a+[927^33 + (ai + n23)7tV - (03 - n^2)'R\ + K^ii - ^22 + n33)7^^] - 

d_[dsn'2 + (ai - 7123)^'- - (02 + nl3)7^^ - i(7lii + 7122 - ^33)^^] " 
[(92 - 02 - ni3)S2 - (^3 - 03 + "'12)-E3 + V?>niiE+ + (7l22 " ?T'33)-B_ + 

2n23£^l]7^^ + [(9i - 2ai + 27123)^2 - 293£;_ + (03 + 7112)^- - 

(3a2 + 57ii3)£;i - |(7iii + 37122 " nzz)E^\R}2 + 

[(9i + 2ai + 27123)^3 - 292^;- - (02 - 7113)^+ + (3a3 - 57ii2)£;i - 

\{nil - 7122 + 37^33)£^2]7^^3} , (A.2) 

^^2 = ^ {[2(9(1 - 3a(i)c72) + (93 - 3a3)CT_ - ni2(T+ - (7111 - 7122)0-3 + ^230-2- 

7ii3(Ti]7^^i - 2[(9(2 - 30(2)0-3) - :^(9i - 3ai)o-+ + 7123^-- - \{n22 - 7^33)0-1 + 
1(7113(73 - 7^l2(72)]7^^ - ^-[2(9(i - 3a(l)7^l3) + (^2 - 3a2 - n^z)'R}- + 
(nil - n^z)n^2 + r^l27^'l - n2zn\] - ^/2,E+[{d^ - 02 - 37ii3)7e^- + 

(93 - 03 + 37li2)7^\ - (9i - ai - 37l23)'^^3 - (^22 " 27lii - 27^33)7^S] + 

^[(92 - 3a2 + 37^l3)a+]7^3_ + [a ^ E ^ U^)} + 

I {(92 - as - 3ni3)(V3Q+ - P-) - (93 - 03 + 87112)^1 + (9i - oi - 3n23)7'3- 
(7222 - 27in - 27133)^2} - \ {3(92^+)7ei_ + 9^+7^^ + lla+TZS} + 
^ {a_[937^^ - 2(a2 + nig)?^^. - (ai - 7^23)7^'\ + i(7Zii + 7122 - n^z)n\] + 
d+[d2n^- - 2(03 - 7ii2)7e3i + 2(ai + 7123)^^ + 2(7iii - 7122 + nzi)n^2] + 
^(7+[(39i - ai - 37^23)7^^3 - (^2 + 02 - 57^l3)7^=^- - (93 + 03 + 57^l2)7^^ - 
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(nn + 3r^33)7^'2] - [^3^- - (^2 + ^a^ + 2^3)^1 - 2(03 - riia)^- - 
(3ai - 57223)^2 + i(3riii + ^22 - ^^33)^3]7^^ + 

[{di + ai - ^23)-^! - (^3 - as - ?^i2)-E3 + 2^13^2 + (nn - ^22 - ?^33)-E- - 
U^n + - r^33)^+]7^^2 - [diE- - {d^ - 2a2 + 2^13)^3 + {a^ + n23)E+ + 
(3a3 + 5ni2)^2 - - ^22 - 3nss)Ei]n\} , (A.3) 



^^3 = - ^ {[2(^(1 - 3a(i)a3) + {82 - 3a2)a+ + ni^cr^ + {nn - n3s)a2 + nucri- 
n23crs]n^i - 2[{d(2 - 3a(2)cT3) - 75(^1 - 3ai)o-+ + n230-_ - |(?222 - ^33)0^1 + 
i(ni3a3 - ni2a2)]n^2 - E_ [2(9(1 - 3a(l)7^l2) - (^3 - 3a3 + r^l2)7^l_ - 
(nn - n22)n\ + n2zn\ - unTZW + VlE+[{ds - 03 + 3nl2)7^'- - 
{82 - a2- 3nl3)7^^ + (^i - ai + 3^23)7^^2 - (?^33 - 2nii - 27222)7^^3] - 
^[(^3 - 3a3 - 3r^l2)^+]7^=^- + {a ^ E ^ TZ^)} - 

1 {(^3 - as + 3ni2)(v^Q+ + P-) - {82 -a2- 3nis)Vi + {8^ - + 3^23)^2- 
(^33 - 2nu - 2^22)^3} - 1 {3(a3^+)7^l_ + dE+nS + lla+n';} - 
^ {a482n\ + 2(a3 - ni2)7t3_ - (ai + r^23)7^=^2 - K^n - ^22 + ns3)n%] + 
a.[8sn^. + 2{a2 + ni^)n\ - 2(ai - n23)n\ + 2(nn + n22 - n33)7^='3] + 
^^+[(3^1 - ai + ?>n2z)V?2 + (^3 + as + 5r^l2)7^'■^- - {82 + aa - 5r^l3)7^^ + 
(nn + 3r^22)7^'s] - [^2^+ + (^3 - 2a3 + 2ni2)^i + 2(a2 + ni^)E_ - 
(3ai + 5^23)^3 - ^(3nn - ^22 + n^'i)E2\R}i - 
[81E+ + (^3 + 2a3 + 2ni2)^2 + (ai - ^23)^- + (3a2 - 5^13)^3 + 
i(nn - 3^22 - n33)Ei]7?.^2 + [(^i - ai - 7223)-^! - {82 - a2 + ni3)E2 + 
2ni2^3 - (nn - ^^22 - riss)^- - |(^ii - ^^22 + ^^33)^-]7^S} • (A.4) 
Expressions ( |A.1[ - |A.4| ) show explicitly that J-'^ab = when 71' = . 
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